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Abstract 

A systematic procedure for performing holographic renormalization, which makes 
use of the Hamilton-Jacobi method, is proposed and applied to a bulk theory of 
gravity interacting with a scalar field and a f7(l) gauge field in the Stiickelberg for- 
malism. We describe how the power divergences are obtained as solutions of a set 
of "descent equations" stemming from the radial Hamiltonian constraint of the the- 
ory. In addition, we isolate the logarithmic divergences, which are closely related to 
anomalies. The method allows to determine also the exact one-point functions of the 
dual field theory. Using the other Hamiltonian constraints of the bulk theory, we de- 
rive the Ward identities for diffeomorphisms and gauge invariance. In particular, we 
demonstrate the breaking of U (1)_r current conservation, recovering the holographic 
chiral anomaly recently discussed in hep-th/0112119 and hep-th/0202056. 
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1 Introduction and Summary 



During the past four years, we have learned much from the AdS/CFT corre- 
spondence [1,2,3] about conformal field theories (CFTs). The basic notion of 
the AdS/CFT correspondence is that a gravitational theory (such as String 
Theory) living on (d + l)-dimensional anti-de Sitter (AdS) space (the bulk 
space) is dual to a CFT living on its conformal boundary. More generally, 
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asymptotically AdS domain wall solutions of (d + l)-dimensional gravity cou- 
pled to certain matter fields are the duals of deformations of CFTs either by 
the addition of relevant operators to the CFT Lagrangian or by the choice of a 
non-conformal ground state involving vev's of certain operators. Using a pre- 
cise recipe one is able to obtain the correlation functions of these (deformed) 
CFTs from the dynamics of the bulk theory. The most important quantity is 
the on-shell action of the bulk theory, which, after suitable regularization and 
renormalization, is identified with the generating functional of the boundary 
field theory. For a recent exposition of the AdS/CFT correspondence including 
an extensive list of references, we refer the reader to the lecture notes [4]. 

The occurrence of divergences in the bulk on-shell action was noted already in 
the earliest AdS/CFT calculations [2,3,5,6]. The first divergence to be explic- 
itly removed from the on-shell action by adding a counterterm was the bound- 
ary volume divergence in the context of pure gravity on an AdS background 
[7]. Soon after, and still for the case of pure gravity on an AdS background, 
the general structure of the divergent terms and the relation of the logarithmic 
divergence to the conformal anomaly of the boundary CFT was discussed in 
[8]. The addition of counterterms also provided a way to give a meaning to 
the notion of energy in asymptotically AdS spaces in terms of the renormal- 
ized Brown- York stress energy tensor [9] without the need of reference spaces 
[10,11,12] 1 

A systematic development of holographic renormalization for bulk gravity cou- 
pled to scalar fields was first given in [15]. This method, which we shall refer 
to as the standard method of holographic renormalization, involves the can- 
cellation of all cut-off related divergences from the bulk on-shell action by the 
addition of counterterms on a cut-off boundary hypersurface and the subse- 
quent removal of the cut-off. The counterterms removing power divergences 
are fully covariant expressions of fields living at the cut-off boundary, whereas 
logarithmic counterterms depend also explicitly on the cut-off. This dependene 
breaks some bulk diffeomorphisms and yields, e.g., the trace anomaly [8]. Most 
recently, the standard method was summarized very clearly by Bianchi, Freed- 
man and Skenderis (BFS) [16] (see also the lecture notes by Skenderis [17]) 
and applied to domain wall bulk geometries dual to deformed CFTs, where 
also bulk gauge fields were considered. Holographic renormalization has been 
used in the calculation of two-point functions in deformed CFTs [18]. 

In a parallel development, de Boer, Verlinde and Verlinde (dBVV) [19] advo- 
cated the Hamilton- Jacobi (HJ) approach in order to separate terms in the 
bulk on-shell action, which can be written as local functionals of the bound- 
ary data. The remaining, presumably non-local, expression was identified with 



For example, compare the treatments of the AdS-Schwarzschild black hole using 
pure AdS as reference space [13] and using counterterms [14]. 
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the generating functional of a boundary field theory. This approach does not 
correspond to the standard one, in particular, because dBVV's boundary field 
theory lives on the cut-off boundary, and because the generating functional 
still contains logarithmic divergences. Nevertheless, dBVV's method remains 
intriguing for its simplicity — it yields formally the correct gravitational anoma- 
lies and provides a remarkably simple bulk description of the renormalization 
group flow in deformed CFTs. A first attempt to use the Hamilton-Jacobi 
method for performing holographic renormalization appeared in [20], where a 
list of references to further studies of dBVV's method can be found. 

Despite its simplicity, the use of dBVV's method to perform holographic renor- 
malization does no appear to be very popular. The main drawbacks in its orig- 
inal formulation seem to be the non-uniqueness of the solutions for the local 
terms to be used as counterterms and the failure to obtain the logarithmic 
counterterms. It is our intention in this paper to overcome these difficulties by 
explaining exactly how the ambiguities can be removed from dBVV's method 
and how logarithmic counterterms are obtained. The main steps of our ap- 
proach, which we shall call the HJ method of holographic renormalization, 
will be as follows. In the first step, we shall use dBVV's method to determine 
the counterterms for the power divergences. Ambiguities in the descent equa- 
tions are removed at the top level by comparison with simple and known free 
field calculations. In the second step, we continue to analyze the hamiltonian 
constraint with two important results. First, we prove that all power diver- 
gences have been removed by the first step, and, second, we obtain explicitly 
the logarithmic divergence and remove it. Hence, we prove that the HJ method 
qualifies as a consistent procedure for performing holographic renormalization. 
We point out that the difference between the standard method and the HJ 
method does not stem from using different renormalization schemes, but re- 
gards only the procedure in which the (exactly identical) counterterms are 
obtained. In both methods a choice of renormalization scheme is reflected in 
the possibility of adding finite counterterms. 

A very interesting aspect of the HJ formalism is the way holographic Ward 
identities emerge. In fact, local symmetries of the bulk theory are encoded 
in hamiltonian constraints for the boundary data. Exploiting the fact that 
the on-shell action is identified with the field theory generating functional, 
these can be naturally interpreted as Ward identities [21,22]. The hamiltonian 
constraint associated with radial evolution can be seen as the Ward identity 
for conformal symmetry, whose anomalous contributions are related to the 
logarithmic divergences of the on-shell action [8]. Similarly, we will see that 
logarithmic terms are also responsible for an anomalous term in the current 
Ward identity. 

Let us give an outline of the rest of the paper and summarize our results. In 
Sec. 2 we shall briefly review the standard method of holographic renormal- 
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ization in order to familiarize the reader with the issue and the complexity of 
the method, and to allow a direct comparison with the approach that we will 
propose. We will then turn to the HJ method with the purpose of performing 
the holographic renormalization. The method will be compared to the stan- 
dard method at appropriate points throughout the paper. For simplicity, in the 
first part of the paper we shall confine our attention to a bulk gravity theory 
coupled to a single scalar field. In Sec. 3 we focus on the procedure of regu- 
larization and renormalization of the bulk on-shell action. This is carried out 
solving the hamiltonian constraint for radial evolution in a recursive fashion. 
Following dBVV, we start with a general ansatz for the leading local terms 
of the on-shell action, which should be covariant and contain all power diver- 
gences. The constraint will then fix all of them, provided some little input is 
given for the lowest levels (Sec. 3.1). As was already discussed in [19,20], there 
might be some obstruction to the solution of the constraint, which is closely 
related to the occurrence of a conformal anomaly. We will make this observa- 
tion more precise and relate the anomaly to the logarithmic divergences of the 
on-shell action, which we can explicitly determine and remove. This is carried 
out in Sec. 3.2, where we also prove that all divergences have been removed. 



In Sec. 4.1, we shall discuss how to determine exact one-point functions in our 
method. By "exact" we mean one-point functions in the presence of sources 
turned on, thus, in principle, encoding all higher point functions of the field 
theory. This is one of the main results of BFS's approach, and any alternative 
method should not fall short of it. Our results will also explicitly demonstrate 
the scheme dependence of local terms in the exact one-point functions. As 
a further application, we consider one-point functions in bulk backgrounds, 
which are called holographic renormalization group flows (Sec. 4.2). The dif- 
ference between operator and vev flows, which describe deformations of the 
dual CFT by either operator insertion or a non-zero vev, respectively, will 
become explicit. Moreover, the most natural choice of counterterms will cor- 
respond to a supersymmetric renormalization scheme. 



In Sec. 5 we will include a U(l) gauge field in our treatment, which should 
correspond to the /^-symmetry current of the dual field theory [23,16]. A 
repetition of the arguments previously explained will allow us to show how 
to solve for the on-shell action up to level four (with some simplifications) and 
to obtain new contributions to the conformal anomaly. We will obtain some 
one-point functions for the new fields and, in particular, we will give a full 
holographic derivation of the complete Ward identities, recovering also the 
chiral anomaly recently discussed in [16], [24] and [25]. 
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2 Holographic Renormalization — Standard Approach 



To start, let us summarize the standard method of holographic renormalization 
as described by BFS [16]. For later convenience, we shall use a bulk metric of 
the form 

ds 2 = dr 2 + gij(x, r)dx l dx^ , (1) 

where i,j = 1 . . .d. The metric (1) is related to the Fefferman-Graham met- 
ric for asymptotically AdS spaces [26] by a simple change of variables. The 
asymptotic region is r — > oo, where the bulk approaches AdS spacetime with 
characteristic length I, i.e., 



The standard method starts by solving the equations of motion in the asymp- 
totically AdS region, which yields the bulk fields in the form of asymptotic 
series like 2 



F(x, r) = e x ^ 1 [fo(x) + e r ' l h{x) + e 2r/l f 2 {x) + ■■■ 
'fo(x)+e r / l f 1 (x)+e 2r / l f 2 (x) + - 



_|_ gA 2 r// 



(3) 



consisting of the two independent asymptotic series solutions of the equations 
of motion. These series expansions are similar to the Fefferman-Graham ex- 
pansion of the metric for asymptotically AdS spaces [26] in pure gravity. 

The assumption that the bulk space be asymptotically AdS implies that the 
fields are asymptotically free, and thus the coefficients Ai and A2 are related 
to the asymptotic mass of the field F, whereas /o and /o are independent 
integration constants. Without loss of generality, we shall assume that Ai > A 2 , 
so that the first line in eqn. (3) represents the leading series solution. Then, f 
is called the source of F. The coefficients /1, f 2 , etc., and similarly, /1, f 2 , etc., 
are obtained by recursive analysis of the bulk field equations in the asymptotic 
region and depend locally on the fo and fo of all fields, respectively. In the dual 
boundary field theory, the source f is the generating current of the operator 
that is dual to the bulk field F, while f is related to the vev of this operator 
(see Sec. (4.1)). 



2 Although this form is not generic, it suffices to describe the standard method. The 
expansion (3) will in general contain not only powers of e r /' and the logarithmic 
terms described later, because bulk interactions and generic values of the confor- 
mal dimensions will generate non-integer exponents. The literature on the standard 
method [15,16] starts with expansions of the form (3), because it covers most inter- 
esting cases, where conformal dimensions are integer or half-integer. However, there 
is no obstruction in principle to apply the method to other cases [27]. 
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For example, for the metric, gij, the coefficients are Ai = 2 and A2 = 2 — d. 
In the case of a scalar field, 0, of bulk mass m, they are Ai = — (d/2 — A) and 
A2 — ~ (d/2 + A), where A = \Jd 2 /4 + m 2 l 2 . The bulk field <fi is the dual of a 
scalar operator with conformal dimension A = d/2 + X. We shall consider only 
matter fields dual to relevant operators, for which Ai < 0, so that they tend 
to zero in the asymptotic region. 3 For even values of d, there are logarithmic 
terms in the leading series for the metric (proportional not to an exponential 
of r, but to r), starting with r x exp[(2 — d)r/l]. Similarly, logarithmic terms 
occur in the leading series of scalar fields for integer A. 

Having obtained the asymptotic form of the bulk fields, one proceeds to cal- 
culate the regulated on-shell action, which is the on-shell action for a bulk 
spacetime with a cut-off boundary at r = p, 



S Teg [fo,p\ = J d d x^T% 



e^W/o) + e^Vt/o) + • • • + ^d(f ) + 0(1) 



(4) 

where the are positive numbers. There are a finite number of divergent terms 
in the limit p — > 00, which comprise the power divergences with coefficients 
cik and a logarithmic divergence with a coefficient a that is related to the 
appearance of logarithmic terms in the expansion (3). The coefficients at and 
a are local functions of the coefficients fo of the leading series of the bulk fields 
and of their derivatives. 

The method of holographic renormalization proceeds now as follows. First, 
the series (3) are inverted for all fields and solved for the sources fo(x) as 
functionals of the F(x,r). Second, the / thus obtained are substituted into 
eqn. (4), which yields the divergent terms of the on-shell action in terms 
of the fields F(x,p) living at the cut-off boundary. Finally, these terms are 
subtracted from S vcg (using counterterms) , and in the remaining expression 
the limit p — > 00 is taken. The finite result (expressed again as a functional 
of the sources / ) is the renormalized on-shell action to be identified with the 
generating functional of the boundary field theory. 

BFS's method is clear and rigorous. The four steps of the procedure — the 
asymptotic expansion of the fields, the recursive determination of the coeffi- 
cients, the inversion of the asymptotic series and the substitution of the in- 
verted series into the on-shell action — form an algorithm that uniquely yields 
all counterterms needed to cancel the p — > 00 divergences of the bulk on-shell 
action. Further finite counterterms could be added and correspond to a cer- 
tain choice of renormalization scheme. One should notice the importance of 
using covariant counterterms, hence the need to invert the asymptotic series 
solutions of the fields. Covariance of the counterterms implies the validity of 
the Ward identity of diffeomorphisms of the boundary field theory. In fact, 



One could also include marginal operators, which have Ai = [17]. 
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the counterterms, when expressed in terms of fo and /o, contain finite terms 
that affect the renormalized on-shell action and do not correspond to a choice 
of renormalization scheme, i.e., they contribute to the physical results. 

It is an interesting question whether there exists a method, which determines 
the counterterms directly at the cut-off boundary without performing the 
asymptotic expansion and its inversion. Such a method is provided by the 
Hamilton- Jacobi approach, which we shall turn to in the next section. 



3 Holographic Renormalization — Hamilton-Jacobi Approach 

As mentioned in the previous sections, we propose to use the HJ method ad- 
vocated for the AdS / CFT correspondence by dBVV in order to determine the 
counterterms. In this section, we shall describe in detail how the counterterms 
are obtained, starting with the power divergences in Sec. 3.1 and discussing the 
logarithmic divergences in Sec. 3.2. The relation between logarithmic diver- 
gences and conformal anomalies will become explicit, a fact that has been used 
by Henningson and Skenderis [8] in the derivation of the conformal anomaly. 

In order to keep the presentation simple and to concentrate on the main steps 
of the method, we consider in this section the action of gravity coupled to a 
single scalar field 4 

+ l -jd d x^H. (5) 

The necessary formulae for the HJ method are summarized in the appendix, 
and the vector sector should be omitted here. We shall include the gauge fields 
in Sec. 5. 

Before plunging into the details, let us briefly compare the HJ approach with 
the standard method of holographic renormalization. In the HJ approach, the 
analysis is carried out directly at the cut-off boundary, so that an asymptotic 
expansion and its inversion are not necessary, except for the leading order term 
determined by the behaviour of the free fields. The recursive determination 
of the series coefficients from the field equations is replaced by the recursive 
solution of a set of descent equations. In order to determine where to stop 
(if the system does not break down by itself), a simple power counting using 
the leading behaviour of the fields is sufficient. The descent equations are 
derived from one equation only, which is the hamiltonian constraint, using the 
independence of the boundary conditions in order to split the constraint into 

4 Details of our notation are given in the appendix. For the more general case of 
several scalars, see [20]. 



x\ g 
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functionally independent terms. One should expect ambiguities, because this 
equation does not contain all information about the field dynamics. In fact, the 
descent equations do not yield a unique solution for the counterterms, but, as 
we shall see, the ambiguities appear right at the top of the descent equations 
and can be easily removed with little input from the (known) results for free 
fields. The subsequent analysis of the descent equations yields a unique answer 
for the divergent counterterms. As a very nice feature, the other constraint 
equations translate into Ward identities of the dual field theory. 



3. 1 Power Divergences 



In this section, we shall briefly outline how the power divergences of the on- 
shell action are obtained using the method of dBVV. In order to separate the 
power divergences, the on-shell action is split up as 

s = s n + s [2] + ■■■ + s [2n] + r , (6) 

where we have denoted power divergent terms by S^k] , k — . . . n. According 
to eqn. (A. 21), the momentum n naturally splits into 

= 7T[ ] + 7T[ 2 ] H h 7T[ 2n ] + 7T r , (7) 

and similarly for g*. We shall make the following premises regarding the coun- 
terterms S[2k]- 

(1) The counterterms are covariant (and gauge invariant, if vector fields are 
involved) local expressions in terms of the metric gij and the field <fi at 
the cut-off boundary as well as their derivatives. 

(2) The term 5W] contains exactly k inverse metrics. (This is an arbitrary 
assumption, but turns out to be very useful for the bookkeeping.) 

(3) The counterterms should completely contain the power divergences of S. 

(4) The counterterms must be universal, i.e., they must contain the power di- 
vergences of S for any asymptotically AdS solution of the bulk equations 
of motion. 



The premises 1 and 2 imply that we can write 
S [0 ] = Jd d x^U(<P), 



Si 



[2] 



M(<i>)g ij d i (l>d j (l> + $(<i>)R 



, etc. 



(8) 
(9) 



The number n of power divergent terms is dictated by premise 3: By power 
counting using the leading asymptotic behaviour of the fields for r — > oo 
[cf. eqn. (3)], one can determine which covariant boundary integrals are gener- 
ically divergent. All of these have to be included. For example, an Spk] without 
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derivatives of matter fields behaves generically like exp[(d — 2k) p/ 1] from the 
metric contributions. Thus, we must include all of these with k < d/2. In 
eqn. (6), T may contain a logarithmic divergence, but shall be regarded as 
finite for the power counting. Finite local terms can be arbitrarily shifted be- 
tween T and the counterterms. This reflects the usual ambiguity of choosing 
a renormalization scheme. 

The main point of dBVV's method is to analyze the Hamiltonian constraint 
[eqn. (A. 17)] by splitting it into a derivative expansion, 5 

H = H[o] + H[2] + h H[4n] + 7~Lr = , (10) 

where H\2k\ denotes those terms in 7i that stem only from the counterterms 
and contain a total of k inverse metrics. Notice that eqn. (10) only re- writes 
the hamiltonian 7i in a form, which makes the contributions from the various 
counterterms explicit. By construction, a certain number of Tipk] are asymp- 
totically stronger than Ti,r, and these can be found by simple power counting. 
One should now try to solve the constraint separately for these (stronger) 
Tt[2k], starting with Tl[o], which leads to a system of descent equations. This 
approach is justified by the universality premise. 

Let us demonstrate the procedure by considering the level and level 2 terms 
of the Hamiltonian. For more details, see [19,20]. The level descent equation 
is given by 

n [0] = -^ju 2 + \ (u'f -v = o. (ii) 

This equation is regarded as a functional equation for U((f>), where V((f>) is 
given. Solving eqn. (11) is impossible in closed form for a generic potential 
V((p) (see [30,31] for a numerical analysis). However, we are interested in 
removing the p — > oo divergences, and thus an analysis in the asymptotic 
region is sufficient. In fact, because the leading behaviour of a relevant scalar 
is (p(x,r) ~ <fio(x)e Xir with Ai < [cf. eqn. (3)], for any finite, but other- 
wise arbitrary value of the source 0o, we can choose the cut-off p sufficiently 
large such as to bring the Dirichlet value <j)(x, p) within an arbitrarily small 
neighbourhood of the fixed point, = 0. The same is true for any number of 
relevant matter fields. Hence, it is sufficient to solve 7Yr ] = close the fixed 
point. For this purpose, we expand it into a Taylor series. 6 Thus, we write 



Such kinds of expansions were proposed previously in the literature as a method 
for solving the Hamilton-Jacobi equations of General Relativity in the presence of 
some matter fields [28,29]. 

6 Note that this is not equivalent to a large p expansion of the descent equations, 
because it does not take into account the different large p behaviours of different 
scalars [27]. 
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V(<f>) as 

V = -^P 1 + \™ 2 <P 2 + ^ 3 3 + + O(0 5 ) , (12) 

where the constant part represents the negative cosmological constant, and m 
is the bulk mass parameter for 0. Of course, this assumes that the potential V 
can be expanded in a Taylor series, which we shall assume to be true. Similarly, 
we expand U as 

U = U + Ul (f)+ l -U^ 2 + ^M 3 3 + ^«40 4 + O(0 5 ) . (13) 

Now, we obtain from eqn. (11) a set of coupled equations, from which the 
coefficients Uk should be recursively determined. It turns out that there is no 
unique solution, but input from the (known) free field behaviour will fix the 
ambiguities. First, we see that the action (5) does not contain terms, which 
are linear in 0. Thus, there cannot be divergent terms linear in 0, and we shall 
set 

ui = . (14) 
This solves the term of order of eqn. (11). Then, we find from the constant 
term 

M 2 = . (15) 

Here, we must choose a sign, and in order to choose the correct one, let us com- 
pare with a simple known case. Liu and Tseytlin [7] determined this leading 
counterterm as 

d ~ 1 (, R \ 

which cancels the volume divergence of the on-shell action for pure gravity [7]. 
This is easy to verify for a bulk AdS spacetime with a flat cut-off boundary. 
Next, the quadratic term of eqn. (11) yields 

(u 2 ) 2 + ju 2 - m 2 = , (17) 

where eqns. (14) and (16) have been used. Again, the solution is not unique, 
and we must resort to the free scalar field in AdS background, whose on-shell 
action is 



S = -Jd d x^<pd r <P, (18) 
and whose leading behaviour (obtained from the free equation of motion) is 

<f)(x,r) ^e { - d/2+x)r/l 4>(x) (19) 

with A = \Jd 2 /A + m 2 l 2 . Hence, the leading behaviour of the on-shell action 
(18) is 

5 ~27/^ 02 H + A ) • (20) 
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This tells us that the correct choice amongst the solutions of eqn. (17) is 

"» = tH +a ) ■ (21) 

Continuing with the cubic term of eqn. (11), we obtain 

so that, for A 7^ d/6, we uniquely find 

US = 6A^ 3 ■ (23) 

In contrast, for A = d/6, the coefficient M3 remains undetermined, and we have 
a remaining term, 

Kcm = -^30 3 + O(0 4 ) . (24) 

We shall postpone the discussion of 7Y rcm until Sec. 3.2. The fact that U3 
remains undetermined is a sign of the renormalization scheme dependence. In 
fact, 

d d x^g<f ~ e dr/l [ e (-<V2+AW*j 3 = 1 for A = d/6 , (25) 

so that the undetermined counterterm is finite. More generally, a breakdown 
of the level descent equation will occur, if 

for some integer k > 2 [31]. In this case, the coefficient Uk remains undeter- 
mined, and the boundary integral of <p k is finite in the p — > 00 limit. A relation 
similar to eqn. (26) can be found also at the higher levels. 

For later use, let us also consider the 4 term of the constraint (11), which 
becomes 



241^" ~' 24 4(d-l)P\2 "J ' 2(6X-d) 2 



-±(4A-d)--i-— -— - --A +—-^— = (27) 



after the coefficients u , U\, u 2 and w 3 have been inserted. Clearly, if A 7^ d/4, 
eqn. (27) determines w 4 uniquely, but, if A = d/4, the descent equation breaks 
and there is a remainder 

*-= (-s- 645^ + 2 S i )* 4+ °^- (28) 

It turns out that the relevant contribution to l~C renL might be zero in particular 
cases. An example is the potential responsible for the GPPZ flow [32], where 
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d = 4 and A = d/4 = 1. The breakdown occurs for k — 4, but we have 
i>4 = — 8// 2 and t> 3 = 0, so that there is no 4 contribution to H rcm . 

It is useful to pause for a moment and reflect on some details of what we have 
done so far. The idea for finding the power divergent terms is to write down 
the most general set of covariant local counterterms up to the necessary level 
and solve for them by asking that they satisfy the hamiltonian constraint for 
arbitrary and independent values of boundary conditions (Dirichlet values at 
the cut-off boundary) of the fields. One might object that this assumes that all 
power divergences can be cancelled by covariant local counterterms, whereas, 
in the standard method, this is shown by explicit construction. At this stage 
of the procedure, we do not have any proof that the procedure applied so far 
yields all power divergences, but this will be proven in Sec. 3.2. 

The origin of the ambiguities we have met becomes clear when considering 
the hamilton equations of motion, e.g., eqns. (A. 13) and (A. 14). The choice 
of sign of «o reflects the fact that one can choose both, r — > oo and r — > — oo, 
as the asymptotic region, and the choice in eqn. (16) corresponds to r — > oo. 
Similarly, the two possible solutions of eqn. (17) describe, via eqn. (A. 14), the 
two independent asymptotic behaviours of the scalar field. The counterterm 
must describe the leading one, because it is asymptotically strongest, and this 
is ensured by the choice in eqn. (21). It is now also clear that further ambi- 
guities can appear neither for higher powers of nor at the higher levels. In 
fact, the hamilton equations of motion provide a one-to-one correspondence 
between the first few terms of the leading asymptotic series and the coun- 
terterms. Hence, since the leading asymptotic series solution is unique for any 
given sources, there is exactly one set of universal covariant counterterms. 

Let us now turn to the level 2 equation, which reads 



By the universality premise, functionally independent terms must vanish sep- 
arately, and this yields the following three descent equations, 




(29) 



d-1 
d-1 




(30) 
(31) 
(32) 



d-1 4 ' 
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These have been analyzed in [20], and the solutions are 

1 



with 



and 



$ = $ O + $ 1 + 1$ 2 2 + O(0 3 ) , (33) 



$0 = -4(^2) f ° r ^ 2 ' (34) 

$i = 0, (35) 

M=^-^ + O(0) forA^l. (37) 

For d = 2 and A = 1 the coefficients <3> , $2 and M remain undetermined, 
respectively, leaving a remainder 

. -iVW^-^ + O^fi 2 ) (A=l), 

!; W + O(0,i? 2 ) (d = 2). { 



Similarly, an analysis of the level 4 descent equations has been carried out 
in [20]. A breakdown at level 4 occurs, if d = 4 or A = 2. For d = 4, the 
unresolved remainder is 

U Tcm = - i^') + 0(i? 3 ) . (39) 



The method can be extended recursively to any desired level, of course with a 
rapidly increasing amount of complexity. The highest level necessary to cancel 
the power divergent terms in S is found by power counting. As discussed 
above, the method yields all covariant boundary terms in the on-shell action, 
which are divergent by power counting. These are interpreted as counterterms 
to be added, although it has not yet been shown that all power divergences of 
S are obtained. It also gives information on possible scheme dependent finite 
terms, and one may choose to add also irrelevant counterterms. We shall turn 
to the analysis of V in the next section, where it is shown that it is either finite 
or logarithmically divergent, which is indistinguishable from finite behaviour 
in the power counting. 



3.2 Logarithmic Divergences and Anomalies 



We have seen in Sec. 3.1 that a zero coefficient in front of one of the unknowns 
of the descent equations leads to a breakdown of the recursion and generically 
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leaves behind an unresolved remainder, 7i rem - We shall discuss the meaning of 
7i rem now. In particular, we will show that it leads to a logarithmic divergence 
of T and to a conformal anomaly in the boundary field theory. If we find that 
H Icm = despite the breakdown, as it happens at level in the case of the 
GPPZ flow potential, then there is no contribution to the anomaly. 

The argument continues with the analysis of the Hamiltonian constraint, 
which, by solving the descent equations up to the breakdown, has been re- 
duced to 

n rcm + n r = o. (40) 

Moreover, we can express Hr to leading order as 

8 

Hr = 8q^q[o]ij - ^— j-<?r<?[o] + ^r^[o] H 

-q r + U'ttt + ■■■ 



d 

1 \2 gij 5T ( d \ ,1 5T 
_ _m3_ 1_ 1_ ^ ^ 

1 [ y/9 $9ij V 2 / Vd °<t> 



+ •••• (41) 



Here, we have used the results for U from the previous section, and the ellipses 
stand for terms that are irrelevant in the p — > oo limit. Defining (remember 
A = d/2 + A) 

2qa 5T , . , 1 5T , A . 

A=^— + (A-d) ( f>—— , (42) 

we find from eqns. (40) and (41) 

A = -lH rem + ■■■ . (43) 

The quantity A is what dBVV call the conformal anomaly of the boundary 
field theory. This identification would follow from eqn. (42), if T were the 
generating functional of the boundary field theory. Formally, A has the form of 
an anomaly, as can be seen in the particular cases. However, this interpretation 
is not correct. First, T contains logarithmic divergences and is therefore not 
the generating functional of the boundary field theory. Second, A scales like 
e -dp/i f or p ^> oo, i.e., it has no finite limit, and, third, the boundary data 
gij and are not the (renormalized) sources of the boundary field theory. 
We shall derive below that the true anomaly is obtained from A by a simple 
rescaling and taking the p — > oo limit. However, it is important to notice 
that the boundary integral of A is finite for p — > oo and has the form of an 
integrated anomaly, i.e., it is conformally invariant. 

But first, we shall evaluate the logarithmic divergence of T. The argument uses 
the diffeomorphism invariance of T and is similar to the arguments used in the 
derivation of the conformal anomaly by Henningson and Skenderis [8]. T is a 
functional of the boundary data and and depends explicitly on the value 
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of the cut-off, p. We shall write it as T[gij(x, p), (f>(x, p); p}. By construction, T 
is invariant under bulk diffeomorphisms, 

r[ gij (x,p),<f>(x,py,p} = Y[g' l0 {x\p')^\x\p') ]P '} . (44) 

In particular, under the following change of variables, 

x' = x , r' = r(l + a) , (45) 



the variation of V becomes to first order 



T[gij(x,p'),(f)(x,p');p'] - V\gij{x,p),(f>{x,p);p] 

\r AFl 

+ 0(a 2 ) . (46) 



5T 5T 



= pa J d d x 

Hence, after expressing d p gij and d p <j) to leading order, we find 

d p T = - l Jd d x^A+-- - , (47) 

where the boundary values and <fi are kept fixed when differentiating with 
respect to p on the left hand side. Since the integral on the right hand side of 
eqn. (47) is finite for p — > oo, this equation implies 

r = j J d d x^gA + finite terms , (48) 

which is the logarithmic divergence of T. Eqn. (48) proves that T does not con- 
tain power divergent terms. In other words, we have found all power divergent 
terms of the on-shell action by the analysis in Sec. 3.1. 

It is interesting to ask why a breakdown of the descent equations is related 
to a logarithmic divergence. In the standard method, one solves the equations 
of motion in form of an asymptotic series. It is well known that, for certain 
powers of the leading behaviour, the series form of the solution will "break 
down" , and one requires logarithms in order to continue the recursion. In the 
HJ method, the descent equations determine the momenta as a series of terms 
that locally depend on the fields, which, by virtue of the hamilton equations 
of motion, corresponds to the leading asymptotic series solution. Hence, the 
descent equations break exactly at the point, where the asymptotic series in 
the standard method would require a logarithm, and this translates into a 
logarithmically divergent term in the on-shell action. 

Having found the logarithmic divergence, we define the renormalized on-shell 
action as 

S rcn = Urn (r + pj d d x^gH rcm ^ . (49) 
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Since A has the form of a conformal anomaly, but vanishes for p 
natural to define a finite, rescaled quantity (cf. [20]), 



A = lim (e dp/l A) . 



oo, it is 



(50) 



A is formally identical to A with and <p replaced by g^ and (f> (the coeffi- 
cients of the leading terms in the asymptotic expansions), respectively Then, 
the expression 

J d d x/§A = Jim J d d x^A (51) 



is conformally invariant, and we confirm that A is the conformal anomaly of 
the boundary field theory, 7 



<T)-(A-d)*<0> = 



V0 Vfi 1 50 



= lim 

p— »oo 



^ * + (A-d)*-L-i' 



lim fe^U) = i. 

p^oo \ / 



(52) 



In passing from the first to the second line we have used the fact that the 
integrated anomaly is conformally invariant. 



3.3 Special Case A = 



There is one logarithmic divergence that has not been addressed so far and 
which occurs in the case A = 0. In this case, no breakdown of the level 
descent equation occurs, although the boundary integral of 2 is finite. Hence, 
this divergence must be derived in a different fashion. This was considered 
also in [20], but the argument there is not complete. 

In the case A = 0, the asymptotic behaviour of the scalar field is not given by 
the generic expression (3), but follows 



Thus, 



4>(x, r) 



8 r 



-dr/(2l) 



~4>{x) + <j)(x) + 



(53) 



(54) 



7 We take —<j) as the source coupling to the operator O. The minus sign differs from 
the usual picture in the literature, but emerges from an improved correspondence 
formula [33]. Moreover, our energy momentum tensor differs by a minus sign from 
the convention of BFS. 
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but we have also 

d r (f) = 7T = U' + 7r r + • • • . (55) 
Comparing eqns. (54) and (55) and using eqns. (13) and (21), we can read off 

7r r = + , (56) 

where the ellipses stand for terms that are exponentially suppressed compared 
to the one written. In order to find the logarithmic divergence in terms of the 
fields living at the cut-off, write 

= 7r r ~ -0 , 57) 

from which follows that 

T~±-Jd d x^g<l>\ (58) 

in addition to any logarithmically divergent terms discussed in Sec. 3.2, which 
might arise for even d. Hence, in the case A = 0, the renormalized on-shell 
action is defined as 



5 ren = Urn (t + pj d d x^H rcm ~Y p J d d x^j<f) 2 



(59) 



The other special feature of this case is the calculation of the anomaly. As 
in Sec. 3.2, we can use the fact that, by solving the descent equations, the 
hamiltonian constraint has been reduced to eqn. (40), but Hr now has the 
form 

1 



■y?r9[o] + vr r 7r [0 ] + -vr r + 



(60) 



We have to include the term 7Tp, because in the power counting it has the same 
strength as the other terms. This yields 



(T) + d -HO) 



S r , 



2§ij 5 1 5 

y/g 5gij 2 

29ij__S d _L 5 




lim 

p— »oo 



(61) 



As before, we have used the conformal invariance of all logarithmic countert- 
erms. In eqn. (61), A is defined by eqn. (43) [eqn. (42) does not hold here], 
so that A denotes the gravitational conformal anomaly for even d, and the 2 
term is the matter conformal anomaly. 
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4 One-point Functions 



4-1 Exact One-point Functions 



One of the main results of the standard method of holographic renormalization 
is that it yields formal expressions for exact one-point functions [16]. Exact 
one-point functions depend on finite sources and thus contain information 
about all higher point functions of the theory. It is therefore desirable to obtain 
the same results using the Hamilton- Jacobi approach. This is straightforward, 
given the renormalized on-shell action, eqn. (49) [or eqn. (59) in the case 
A = 0], and we shall demonstrate the procedure by considering the scalar 
one-point function for the cases d — 4, A = 1 and A = 0, d arbitrary. 

Generally, we obtain 

1 ^iSVnn 



(0(x)) 



Vg 5<j>(x) 



= — lim 

p— >oo 

= — lim 



(d/2+X)p/l 1 



^fg5(f)(x,p) 
SH r 



T + pj d d x^H t 



e W2+x) P /ii r + p 



8(f) 



(62) 



but, in practise, it is more useful to substitute 7r r from the split form (7) of 
the hamilton equation of motion (A. 14), 

7T r = d r (f) - 7T[ ] 7T[ 2n ] , (63) 

because at this stage 7i~r is not known explicitly in terms of the source. More- 
over, it is necessary to write down an asymptotic expansion for beyond 
leading order, namely of the form of eqn. (3). However, in contrast to the 
standard method, we do not have to solve for the subleading coefficients using 
the equations of motion, an ansatz with the subleading /1, f'2, etc. undeter- 
mined is sufficient. In practise, finding a general ansatz asks for some care, 
because generic bulk interactions might generate all sorts of exponents in the 
subleading terms. However, for our purposes, we shall regard this step as done. 

Specializing to d — 4 and A = 1, the potential V has the form 

This example is a generalization of the potential generating the GPPZ flow 
[32]. According to the results of Sec. (3.1), the function U(4>) occuring in the 
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leading counterterm has the form 

w)=4-^ + > + > 4+ -- (65) 

Remember that the coefficient remains undetermined. Similarly, we have 
for the functions $(0) and M(0), 

= ~l + \® 2<p2 + " ' ' (66) 
M(0) = m + • • • , (67) 

with $ 2 and mo undetermined. The relevant asymptotic expansion for is of 
the form 

r ) = e-^0 + e- 2r /Vi + e" 3r/ V 2 + ye" 3 ^ + • • • , (68) 
while it is sufficient to use the leading term for the metric, 

g i3 (x,r)=e 2r / l g lJ + -- - . (69) 



The unresolved remainder from the descent equations is given by eqns. (28), 
(38) and (39), 



-T7 ( ^ - R *i Rl3 ) - ^V>V,0 - i-R(f) 2 + 



16 V3 



'fv! 



12 



1 1>4 

3/2 ~ 24 



Thus, after substituting everything into eqn. (62), one finds 



(70) 



(°>='™{ eP/ '(T 01 + T* 2 ) + 



2 , 



v 4 



/ 2 ?; 2 \ « 1 

3 I 53 



6 



6 3/ 2 



+ 702 - T^2 + W01 + 7^0 3 - ™ O V 2 + *2^0 

/ / 6 



(71) 



Notice that the leading divergent term, proportional to e 2p//i , has explicitly 
cancelled. Moreover, since we have removed all divergences from the on-shell 
action (as was proven in Sec. 3.2), the correlation functions must be finite. 
Hence, the remaining two divergent terms in eqn. (71) must vanish by con- 
struction. Thus, we obtain the coefficients 0i and ip2 of the leading asymptotic 
solution, 



01 



l 2 v 3 



2 , 



^2 



I 2 



V 2 - ^0 - 



v 4 



6 + 3/2 



/V\ A 
3 \ 73 



(72) 
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Substituting eqn. (72) into eqn. (71), we finally obtain 



(73) 

where the non-local behaviour is encoded in the undetermined coefficient 02 
of the subleading series solution, which is obtained by imposing regularity of 
the field in the bulk. The arbitrary parameters U4, $2 an d mo represent the 
scheme dependence. In particular, there exists a renormalization scheme, in 
which 02 alone represents the exact scalar one-point function. 



This example underlines the reverse approach of the HJ method with respect 
to the standard method of holographic renormalization. While, in the latter, 
one first determines the sub-leading coefficients in the asymptotic expansion 
of the fields from the equations of motion and then proceeds to renormalize 
the on-shell action, in the former one first renormalizes the on-shell action and 
then obtains the sub-leading coefficients from the finiteness of the correlation 
functions. One could ask whether the subleading coefficients obtained in eqn. 
(72) agree with those calculated in the standard method, i.e., whether they are 
consistent with the equations of motion. In the present example they obviously 
do, as seen by comparison with [16], but the answer is affirmative in general, 
because the coefficients follow uniquely by the procedure given above. 



The case A = is even simpler. Since the renormalized on-shell action is given 
by eqn. (59), we obtain 



(0) = 



J_5_S, 

Vg 54 

lim 

p— »oo 



lim 

p^oo 



dp/{2l) P 



, 1 - -(, 
I V p 



(74) 



In the second line we have used the fact that 7i rem does not depend on in 
the present case, while the last result follows from eqns. (53) and (56). This 
result is in agreement with our expectation, viz., <j> corresponds to the vev in 
the boundary field theory. 



The same procedure can be used to calculate the exact one-point function for 
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the energy momentum tensor, (T^-). The general formula is (for A 7^ 0) 
2 8S ren 



= lim 
= lim 

p— >oo 



s (d-2)p/i 



2 5 



2p 5 



d d x^gH r 



(75) 



For d = 4, the full analysis is quite tiresome and involves also scheme depen- 
dent counterterms, which are quadratic in the boundary curvature. We shall 
not perform the explicit expansion, since, for the GPPZ and the Coulomb 
branch flows, the result is known [16]. We only point out that the trace of the 
energy momentum tensor one-point function is given by the anomaly formulae, 
eqns. (52) and (61), and its divergence satisfies the Ward identity 

Vj(T ij ) + (<9)VV = . (76) 

Eqn. (76) follows from the constraint (A. 18) and the fact that all countert- 
erms are diffeomorphism invariant, because they are expressed as covariant 
boundary integrals. 



4-2 One-point Functions in Holographic RG Flows 



The calculation of one-point functions in holographic RG flows, i.e. those with 
the sources set to their background values, is particularly simple. Holographic 
RG flows are solutions of the bulk field equations of the form 

= 0(r) , 9ij = e 2A (%- , (77) 

which are obtained by solving the first order equations [34,35,36] 

d r( t> = W\4>) , 8 r A = -^L_W{4>) , (78) 

if the potential V can be written in terms of a function W((p) as 

V{4>) = -^jW 2 + \ (W'f . (79) 

In addition, for r — > 00, the solution should approach a fixed point of W with 
negative value, so that the bulk becomes asymptotically AdS. 

It is obvious that eqn. (79) is identical to the level zero constraint, eqn. (11), 
when W is substituted for U . Similarly, eqn. (78) are the Hamilton equations 
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of motion, (A. 13) and (A. 14). This means that W necessarily is a solution 
of the level zero descent equation, but in order to be used in the leading 
counterterm, it must have an expansion 

•"<*> = + s ("5 + A )*" + -- (80) 

In particular, the quadratic coefficient must be —d/2 + A with A > 0. 8 Let us 
assume for the moment that this is the case. Then, after identifying W = U, 
the calculation of the scalar one-point function becomes trivial: 



(O) = - lim 



e W2+x) P /i ( dA _ ir + p 5 J^L 



(81) 



The second equality holds because of eqn. (78) and because neither the higher 
order counterterms nor H Tcm contribute in the Poincare invariant background. 



Similarly, for the energy momentum tensor one finds 



. (82) 



Eqn. (81) explicitly shows that prepotentials W with the expansion (80) are 
flows generated by the addition of a relevant operator to a CFT Lagrangian. 

Notice also that the choice U = W for the leading counterterm means that 
we fix the coefficients in U, which are left undetermined by a possible break- 
down of the level zero descent equation. This corresponds to a supersymmetric 
renormalization scheme, which is supported by the vanishing of the energy 
momentum tensor (82) and of S Icn for a bulk solution of the form (78). 

The second possibility, which we shall now discuss, is that W has the expansion 

^>=-^4(-H 02+ -' (83) 

Here, A is again positive. In this case, the background scalar field has the 
asymptotic behaviour 

0(r) = e -W2+x) P /i^ + . . . ? (g4) 
which is entirely sub-leading with respect to the generic behaviour. Then, 



8 This is the only detail, which is not fixed by the level zero constraint, the asymp- 
totically AdS metric and the existence of the fixed point. 
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eqn. (62) yields 



(0(x)) = - lim 

p—*oo 



— lim 

p— >oo 



e (d/2 + x )P /i ( w _ UIJrp 6 2k 



,(d/2+X)p/l _ 



2A 

T 



+ O(0 2 



2A 

T 



(85) 



Notice that the terms in 7i rem , which do not contain the boundary curvature, 
stem from a breakdown of the level zero equation and involve with at least 
cubic power, and these terms vanish in the p — > oo limit due to the weak 
asymptotic behaviour of 0. Eqn. (85) explicitly demonstrates that holographic 
RG flows with W of the form (83) correspond to deformations of the boundary 
CFT by switching on a vev of the scalar operator. Also in this case the one- 
point function of the energy momentum tensor vanishes. From eqn. (75) one 
finds 



(Tij) = lim {e 



(d-2)p/l 



lim I e dp/l 



p^OO 



{W-U)g l] +2pg l3 0^)\} 

= , 



+ 2 P ^O(0 3 ) 



(86) 



again due to the weak asymptotic behaviour of 



Last, let us comment briefly on the case, where W has the expansion (80) with 
A = 0. The asymptotic behaviour of the RG flow solution is given by eqn. (53) 
with = 0. The scalar one-point function has been calculated in Sec. 4.1 and 
is given by eqn. (74). The energy momentum tensor again vanishes, since, in 
the additional counterterm, the scalar field is not strong enough to compensate 
the 1/p factor. 



5 Contributions of the Vector Sector 



In the following we will add the vector sector of the theory to the gravity and 
scalar sectors discussed so far. We will show how the Ward identities previously 
derived are modified by the presence of interacting vector fields, including ex- 
pected new contributions to the conformal anomaly, and discuss the additional 
current Ward identity associated with the U(l) gauge invariance. In this way 
we provide an alternative derivation of the holographic chiral anomaly, already 
obtained in [16] by using the standard approach to holographic renormaliza- 
tion and in [24] by explicitly analyzing some dual supergravity solutions. 

The physics of an RG flow in the presence of a residual U{1) _R-symmetry is 
encoded in a bulk d+1 dimensional action, which contains a "massive" vector 



23 



field together with a Stiickelberg field a, interacting with scalars and grav- 
ity (see [23,16]). Thus, our starting point is the gauge-invariant action (A.l). 
The Hamiltonian constraints associated with this are derived and summarized 
in the Appendix. Before we begin the analysis, let us say two words on the 
issue of the asymptotic scaling of the fields. When we include the new fields, 
the equations of motion become rather intricately coupled, and we have to 
make sure that interactions between the fields do not influence their leading 
asymptotic behaviour. As a simplification, one can consider gravity and vec- 
tors in some fixed scalar background. This is the assumption made by BFS, 
and we borrow their analysis of the leading behaviour of the vector sector 
fields. From [16], we see that 

B i {x,r) = B i {x) + --- (87) 
a(x,r) = a(x) + --- (GPPZ) (88) 

a(x,r) = j&(x)+a(x) + --- (CB) . (89) 

v 

for r — > oo. 



5. 1 Descent Equations and Conformal Anomaly 



As discussed in Sec. 3, in order to analyze the descent equations obtained from 
the constraint Ti, — 0, it is convenient to group the terms of the local part of 
the on-shell action into different levels. The lowest possible counterterms of 
the vector sector are of level two and read 



QV 

^ [2] 



d d xy/g 



iN^BiBj + Ip^W^B 



(90) 



After computing the momenta stemming from the sum of (8), (9), and (90), 
and inserting them into eqn. (A. 17), one can see that the level zero equation 
is unchanged, and one can always solve the level two contributions by setting 
P{4>) = 0. This is justified, since we do not expect counterterms linear in 
In this way the analysis of the gravity-scalar sector remains unchanged with 
respect to Sec. 3.1. The new equation to be solved is 

(U'N' d-2 N 2 M 2 \ 

We proceed in the usual way by further expanding U, N and A4 2 in powers of 
0. Notice that, for the flows we are interested in, K(<f) = 0(1) and M 2 (4>) = 
0{4> 2 ). Using the solution of U obtained in Sec. 3.1, we first get 
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Again we shall compare to pure AdS/CFT results in order to select one of 
the two possibilities in the above equation. Looking for instance in [37] we 
learn that one should have N = for asymptotically massless vector fields. 
Proceeding up to quadratic order in 0, we unambiguously determine 

iVi = (93) 
N 2 = J^M 2 2 . (94) 

A breakdown occurs, if A = 1, thus giving the contribution 

H Icm = ^Mffi&Bi + O(0 3 , B 4 ) (A = 1) . (95) 

Following the general analysis of Sec. 3.2 this yields a logarithmic divergence 
of T and a contribution to the conformal anomaly for scalar operators with 
A = 1. 

Proceeding to the next levels increases the number of invariants very quickly, 
and a complete analysis already at level four would be extremely tedious. 
The gravity-scalar sector has been analysed in [20], and in the following we 
restrict to the vector sector in a fixed scalar background. A possible basis of 
independent invariants for the level four on-shell action is the following 

{B t B l )\ V i B i B j B j , BiBjVW, F l3 F^ , (V^) 2 , 

v, r>, v j ir . v, r>j v n j . it,, ir n j . v n j it,, . ( ' 

It turns out that the resulting level four descent equations, which we do not 
write here, can be consistently solved setting to zero the coefficients of all the 
terms except for F^F % K This leads to the descent equation 

4 G - ]k ) FijF ij = , (97) 




4/ u 4 

where G is the constant part of the coefficient of the F^F 1 ^ counterterm. The 
descent equation (97) breaks for d = 4, leading to a logarithmic divergence of 
T and to a contribution to the conformal anomaly. 

Hence, we have found the following anomaly contributions from the vector 
sector, up to level four: 

A= [-\IM&B& (A = l), 
[llKoFijfr (d = 4). 

They agree with the logarithmic counterterms of [16], and the d = 4 contri- 
bution was determined also in [38,22]. 
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5.2 Conformal Ward Identity 



Let us show how the anomalous conformal Ward identity is affected by the 
presence of currents. First, notice that the level zero terms of the E l and 7r a 
momenta vanish, because the vector part of the on-shell action starts at level 
two [cf. eqn. (90)]. Therefore, 7ir is the same at leading order as in eqn. (41). 
This implies that exactly the same form of conformal Ward identity holds 
here, though with the additional anomaly contributions derived in Sec. 5.1, 
eqn. (98), namely 

(T) - (A - d) (O) = A . (99) 

This might seem slightly odd. However, it is consistent with the scaling dimen- 
sions assigned to the various fields, which are read as usual from the leading 
r dependence of the asymptotic expansions. In particular, recall that we have 
the following behaviours near the boundary (r — > oo) 

0~ 0e- (d - A)r/l A~Ai a~a. (100) 

This implies that the sources A, and a do not transform under Weyl rescaling, 
and eqn. (99) correctly reflects this fact. On the other hand, it is well known 
that a conserved current should have conformal dimension A = d — 1, and 
thus couple to a vector source of weight one. In the context of the AdS/CFT 
correspondence this is realized by using as source the (rescaled) vector field in 
a local Lorentz frame, A a = e*Ai [37]. Hence, we should have used the frame 
fields A a and the vielbeins e" as "coordinates" in the Hamiltonian treatment 
of the system. This would complicate the analysis unnecessarily, but it is easy 
to correct eqn. (99) for this misuse. In fact, setting S ren [gij, Bj\ = S^nfe", B a ], 
a straightforward application of the chain rule gives 

(f a t )=T a t +B t (J a ) . (101) 

After tracing the above expression and inserting it into eqn. (99), we get 

(T) - (A - d) (O) - B a (J a ) = A . (102) 

Anticipating the next section, we can use gauge invariance to rewrite eqn. (102) 
in its standard field theory form as 

(T) - (A - d) j> (O) - a(O a ) - A a (J a ) = A , (103) 

where we have used the freedom to add a total derivative to the anomaly. 
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5.3 Diffeomorphism and Gauge Ward Identities 



In this section we shall derive the Ward identities that follow from the hamil- 
tonian constraints associated with bulk diffeomorphisms and gauge invari- 
ance. From the constraints (A. 19) and (A. 18) and the premises on the power- 
divergent counterterms stated in Sec. 3.1 we get 



respectively. How these relations translate into Ward identities for the finite 
field theory quantities depends, as usual, on whether logarithmically divergent 
terms hidden in T possibly violate the constraints. Let us assume here that all 
logarithmic counterterms are obtained by the breakdown of the descent equa- 
tions as described in Sec. 5.1. These counterterms are, by construction, both 
covariant and gauge invariant. Thus, the relations (104) and (105) directly 
translate into the following Ward identities for gauge and diffeomorphism in- 
variance, respectively, 



The expression given by BFS as the diffeomorphism Ward identity (eqn. (4.19) 
in [16]) is obtained by substituting eqn. (106) into (107). 

5.4 Special Case 

As in Sec. (3.3), one has to be more careful, if the leading term of a field 
is logarithmic, as is the case for a in the Coulomb branch flow with a fixed 
scalar background [cf. eqn. (89)]. In fact the following analysis applies only if 
we consider gravity and vectors in a fixed scalar background (denoted with a 
bar). That is, we are neglecting the back reaction of a possible scalar source 
turned on. This would complicate the situation considerably, as can be seen by 
inserting the leading behaviour of a from eqn. (89) into the action: The source 
a contributes to the asymptotic mass of the scalar <fi violating the assumption 
that the leading order behaviour does not depend on interactions. 

Let us discuss this special case and derive the anomalous Ward identities for 
the CB flow. Using eqns. (89) and (A. 15) and following similar steps as in 
Sec. 3.3, we find 



V^ + 7T Qr = 
2V^ r - V^Trr - F i3 E 3 T - B^ aV = , 



(104) 
(105) 



V,(J l ) + (O a ) = , 
V J (T 4J ) + (O)VS + (J^Fij + (O a )(A + d t a) = . 



(106) 
(107) 




(108) 
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From this follows that 



4p 



(109) 



in addition to the logarithmically divergent terms stemming from 7i rem - The 
renormalized on-shell action is given by 

S vcn = lim (r + pj d d x^H Icm - ^3 J d d x^gp a 2 ^ . (110) 

Notice that, with the scalar background fixed, the additional counterterm in 
eqn. (59) is irrelevant, and we have not included it here. 

The additional counterterm in eqn. (110) is not gauge invariant, which will 
violate the Ward identity (106). In fact, it follows straightforwardly from 
eqns. (104) and (108) that 



(Ill) 



Hence, we obtain 



(112) 



where for the last equality we have substituted the specific values M\ = 12, 
= -1/V6 for the CB flow [16]. 

As an aside, we can also compute the vev of the operator dual to a for the 
CB flow, 



(O a ) = 



1 5S r , 



— lim 

p^oo 



*dp/l 



pJ__S_ 

1 y/g5a 



r 



Mj 
4p 



J d d x^g$ 



2 a 2 



= — lim 

p-^OO 



M 2 



lim 

p^oo 



2p 
1 
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2 (j) 2 a 



M 2 2 4> 2 a 



-a . 



(113) 



The results (106) and (112) were derived in [16] within the standard approach 
to holographic renormalization and represent the holographic realization of 
the chiral symmetry breaking. As also noticed in [39,24], it is clear that the 
spontaneous breaking of gauge symmetry in the CB case is essentially dual to 
a Higgs mechanism in the bulk, with <p 2 providing a mass for Bi because of 
its weak asymptotic scaling. 
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One might worry that the leading logarithmic behaviour will upset the diffeo- 
morphism Ward identity (107), but the following calculation explicitly con- 
firms it also in this case. From eqns. (105), (110) and (108) follows 



e-^VjT] + J-A^c^V) + e- dp/l d^(0) - ^-Mj^a 2 
Ap 2p 

+e- d ^ l F l3 (.P) - (A. + d^^Mj^a = . 
Using the explicit results (112) and (113) this yields 

V 3 T\ + VMO) + F i3 {P) " AiVj{J j )] + ^jMl^d.a (-a - a 



-dp 1 1 



(114) 



3 -dp/J 



21 

V 3 T\ + V4(0) + F i3 (P) - A&jiJt) + di&(O a )] = . (115) 
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A Hamilton-Jacobi Method 



In this appendix, we shall summarize the Hamilton-Jacobi approach for a 
system of Einstein gravity coupled to a scalar field and a vector field treated in 
the Stuckelberg formalism. Although the summary is sufficiently self-contained 
for the purpose of this article, the reader is referred to standard texts, e.g., 
[40] , for a detailed description of the method. 

The Hamiltonian treatment of gravity involves the time-slicing formalism, 
which assumes that the bulk space-time manifold can be globally foliated into 
hypersurfaces specified by a "time" coordinate. 9 Of course, with Euclidean 

9 For other cases, such as Taub-NUT spacetimes (see [41,42] and references therein), 
the method must be extended. 
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signature, there is no distinction between time- and space-like directions, but 
the method can be applied equally. As for notation, we adorn geometric bulk 
quantities with a tilde and leave those belonging to hypersurfaces unadorned. 
Greek indices, fj,,v, run from to d, latin indices, from 1 to d, and the 
index r is often used instead of the index 0. Our conventions for the curvature 
tensor are R jki = c^Fjj + V i km Y m j i - (k <-> /), Rij = R k ikj . 

For the applications of the present paper (see [23,16]), we shall consider the 
following action: 



S = J d d+l xjg [-Ifl + i^<R0 + + l -J d d x^H 

1 - 



+ / d d+1 xJg 



(A.l) 



Here, V((f>) denotes a scalar potential, which has a local minimum (stable fixed 
point) at = 0. The second integral in eqn. (A.l) is the Gibbons-Hawking 
term, where H is the trace of the second fundamental form of the boundary 
hypersurface. This term is included in order to remove second derivatives with 
respect to the boundary normal from the bulk integral. The vector sector 
contains a vector field and the Stiickelberg field a, with = + d^a 
and F^ = d^B v - d u B fl . 



Writing the bulk metric as 



9»v 



riiU 1 + n 2 rij ^ 



n; 



(A.2) 



9ijJ 



where n and n l are called the lapse and shift functions, respectively, the second 
fundamental form of r = const hypersurfaces is given by 



Hij = -— (drQij - ViUj - VjUi) . 



(A.3) 



Using geometric identities, eqn. (A.l) can be identically re-written as 

S = J d d+1 x^n ]^(-R + H)Hl - H 2 ) + \g^d^ + V{<j>) 
+ J d d+1 x^n ^K^F^F^ + l -M 2 {^)B^ 



10 



(A.4) 



The action (A.4) is invariant under bulk diffeomorphisms and gauge transfor- 
mations, SA^ = —d^5a. We shall gauge fix the quantities 

n = 1 , n* = , A r = , (A.5) 



10 



Notice that R denotes the intrinsic curvature of r =const hypersurfaces. 
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so that their corresponding equations of motion will enter as constraints into 
the Hamilton formalism, 



_SS_ 

Sn 
SS_ 

Sn 1 
SS 

SA r 



H = , 
Hi = 0, 
£ = , 



where (in gauge fixed form) 

H = \ (R + H)H{ - H 2 ) + l -{d r <P) 2 + l -M 2 B 2 r + ^KgVFriF, 
- y^hon.o - l -M 2 B t B* - l -KF l3 F^ - V , 
= \^M H - H l) - drWit - M 2 B r B l - K^FijFrk 



and 



(A.6) 
(A.7) 
(A.8) 



Q = M 2 B r + V\KF ri ) . 



(A.9) 



(A.10) 



(A.11) 



The gauge fixed action reads 

1 



S = J d d+1 x^g [-- (R - H)Hl + H 2 ) + -(<9 r 0) 2 + -g ij d^drf + V 

(A.12) 



1 

2" 



-M 2 B 2 + /•;,/•;., + ^M 2 B t B* y A / / 



1 



1 



where Hij = —\d r gij, B r = d r a and F ri = d r Ai. 

This gauge fixed system can be treated in Hamilton language. The conjugate 
momenta are 



q* = 


1 SS 


= I - IF*) , 


(A.13) 


y /g5(d r g ij ) 


7T = 


1 SS 




(A.14) 


y/g8(dr4>) 


7T Q = 


1 SS 


= M 2 <9 r a , 


(A.15) 


y/g8(d r a) 
1 SS 


E L = 


= KgVF rj . 


(A.16) 


^gSidrA) 



Expressed in terms of the momenta, the constraints (A.6), (A.7) and (Ai 
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become 



n = m - ^ + ^ + + ±m0 (A.i7) 

+ \R- lg ij di<i>d j( j> -V- l -M 2 B t B* - l -KF l0 F^ = , 

Hi = 2V 3 ql - ndi<j> - B-i ir a - F^ = , (A. 18) 

g = 7la + y i E i = . (A.19) 

It is easy to realize that H coincides with the canonical hamiltonian density. 

The bulk theory is defined on a bulk spacetime with a boundary at r = p, 
where p is a cut-off parameter. In the Hamilton-Jacobi formalism the mo- 
menta of the theory are obtained from the on-shell action S as a functional of 
prescribed boundary data, gij(x,p), <f>(x,p), a(x,p) and Ai(x,p), from 

9 ^/gbgij ' (A.20) 

* = -7=Xi- (A - 21) 

where the variation is with respect to the boundary data. Therefore equations 
(A. 17 - A.19) become constraints to be satisfied by S. 
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